Classical back reaction of low-frequency cosmic gravitational radiation ON On =3 G. Dautcourt Max-Planck Institut fur Gravitationsphysik, Albert- Einstein- Institut, Haus 5, Am Muhlenberg, D - 14476 Golm, Germany We study in a Brill-Hartle type of approximation the back reaction of a superposition of linear gravitational waves on its own mean gravitational field up to second order in the wave amplitudes. The background field is taken as a spatially flat Einstein-de Sitter geometry. In order to follow inflationary scenarios, the wavelengths are allowed to exceed the temporary Hubble distance. As in optical coherence theory, the wave amplitudes are considered as random variables, which form a homogeneous and isotropic stochastic process, sharing the symmetries of the background metric. A segregation of the field equations into equations for the wave amplitudes and equations for the background field is performed by averaging the field equations and interpreting the averaging process as a stochastic (ensemble) average. The spectral densities satisfy a system of ordinary differential equations. The effective stress-energy tensor for the random gravity waves is calculated in terms of correlation functions and covers subhorizon as well as superhorizon modes, where superhorizon modes give in many cases negative contributions to energy density and pressure. We discuss solutions of the second-order equations including pure gravitational radiation universes. PACS numbers: 04.30.Nk, 98.80.Hw, 98.80.Cq > oo O o I - 1— < I. INTRODUCTION Many classical aspects of gravitational waves are still poorly understood. Most of them are connected with the fact that the full theory of gravitational radiation as following from general relativity is extremely nonlin- ear. Surprises can be expected, if the nonlinear regime becomes more deeply explored. Pure numerical meth- ods acting on the full field equations [jl) are important but cannot be the only key to the nonlinear regime. They should be supplemented by suitable approximation methods, which allow an analytical or semi-analytical approach. Already in 1964 Brill and Hartle B pro- posed a scheme, which takes the back reaction of lin- ear gravitational waves on the background metric into account. Usually, the Brill-Hartle method (see also [pi, B) is considered as a high-frequency approximation for gravitational radiation. For cosmological application the full spectrum of gravitational radiation including low- frequency (superhorizon) modes must be studied, if its generation and propagation through inflationary stages is considered. We discuss in this article a different inter- pretation of the Brill-Hartle approach, which allows to treat also low frequencies. The Brill-Hartle method requires an average over small scale ripples of the geometry. How averages can be for- mulated consistently in general relativity, is a basic and still not completely solved problem [pi. Our assumption is to interprete the perturbations of the geometry as ran- dom functions as in optical coherence theory B or in the Monin-Yaglom approach to statistical fluid mechanics [J7) . The averages of the Brill-Hartle method are then taken as ensemble averages, not as averages over space-time regions. How stochastic averaging and a random geome- try can be reconciled with general relativity in a rigorous way is also an open question. The problem is not touched upon here, we adopt instead a field theoretical approach in the sense of the Monin-Yaglom treatment. Along this line we discuss the response of the expan- sion rate and of the wave amplitudes to the mean grav- itational field produced by the waves. We assume that the stochastic process for the wave amplitudes is homoge- neous and isotropic. It is not necessary to suppose it to be Gaussian, since only two-point correlation functions are involved at the approximation level considered here. The background is kept as simple as possible, a flat 3- geometry with scale factor a(rj) and conformal time w is assumed. The treatment is completely classical, we do not consider back reaction effects from quantum grav- ity ||, semiclassical gravity |9j or an Einstein-Langevin equation |1(|. With methods similar to those discussed here, the subject was treated in |0| and |P^[ . The paper is organized as follows. In section II the segregation of the field equations into equations for the wave amplitudes and equations for the background field is discussed and, following earlier treatments Q, Q, the effective stress-energy tensor for the random gravity waves is calculated in terms of correlation functions for the wave amplitudes. Depending on the wave spectrum, major contributions to this tensor may come from waves with wavelengths exceeding the horizon distance. Sec- tion III is concerned with gauge problems. The rest of the paper treats solutions of the averaged field equations with various assumptions for the spectrum of gravita- tional waves and for the presence of matter fields. The aim is to discuss effects of back reaction in general, we do not consider the origin of the waves and will also not take observational constraints into account in this article. In section IV we consider the evolution of the gravitational wave amplitudes for different frequency regimes in cosmo- logical models, using a number of simplifications. Pure gravitational wave models ("geons") are the subject of Section V. Its discussion is complicated by the necessity to include a non-linear back reaction term in the wave equation. In the final section VI possible improvements of the approach are shortly discussed. II. AVERAGING THE FIELD EQUATIONS. We take an Einstein-de Sitter model as background metric and add tensor perturbations, writing^ (i,k = 1,2,3) ds 2 -a(rj) 2 dri 2 + g ik dx l dx k with 9ik — fl &ik + hik, hu — 0, h ik _ k — 0. (1) (2) It is appropriate to write the field equations to this metric in the (3+1) form, assuming a perfect fluid as matter in the background metric: r(3) +k 2 - K lk K lk = 16TrGp m , (3) 9 Kl K, hill K, 0, -K'fk/a + i4 3) + KK lk - 2K a K k l = 4nGgik(p m — Pm)- (4) (5) The covariant derivatives are denoted by a prime and taken with respect to gi k , and ' = -4-. Indices are moved with the three-metric gi k . R ik and R^ are the three- dimensional Ricci tensor and scalar for g^ The extrinsic curvature Ki k is given by K ik = ~Ya 9 ' ik = ~ a ' Sik ~2a h ' ik ' Up to second order we have g ik = S ik /a 2 - h ik /a 4 + h u h k i/a 6 . (6) (7) With the decompositions (0) and (0) one can write down the field equations explicitly, including all terms up to the second order. Simplification results from the gauge restrictions for hik, which are taken into account. (3) The Ricci tensor R\ k is given by »( 3 ) 1 h Rik ~ ~^2 hlk > 11 T~ T iLlm\" j ik.lm > i^lmAk ^iLkm ^kLim) + -T—lhl m ,ihl m! k + -^—[(hil im hkl,m — hu jm hkm,l)i (8) and the field equations read explicitly 13 1 + —^(hklhkl,mm + -rhkLrnhkLm ~ -zh k i m h km i) a° 4 2 a' a' 2 1 f h k ih' kl + 3—h kl h kl - —^h' kl h' kl = l&irGp m , (9) 4a 6 - ^-Qhkihki,i gh k ihik,i + ■^-ghkih ik<l -^h kl ti kll + —h' kl h k i,i - -^hki,ih' kl = 0, (10) {— + ~^> dik + ~^2 h ih - ~^ h ik + 2 -X«* a' 2 1 , , a! , a! , -2-jhuhki - 2-±h u h kl + -rh u hu + —gh ki hu a" a' 2a~ 5 ' a" (3) — -^ZK^ ik ^lm^lm ^ QSikhi m hi m + i? ife 4nGgik(ftm-Pm)- (11) The Brill-Hartle method usually starts with the assump- tion that the space-time variation of the small ripples h^ (of order h/X, where A is a typical radiation wave- length) is much larger than the variation of the back- ground metric (of order 1/L). Thus terms in the Ricci tensor which are bilinear in h^p (of order vr) are com- parable to terms involving the background metric (of or- der -jtj), if h rs y. This latter condition excludes low- frequency waves A > L, since h must be sufficiently small to allow a second-order approximation. Terms linear in hjiv,p<r are much larger and should therefore vanish sepa- rately, giving rise to the linear wave equation for h^ v . We do not follow this bookkeeping ( it was criticized in Jig]). Instead, the functions hik are interpreted as random functions. Performing a stochastic average of ([|[nj[ll]) removes terms linear in h ik , but keeps terms bilinear in hik and in the derivatives of hik ■ This allows a segregation of the field equations without restrictions for the wavelengths. The result of the stochastic average can formally be written as (Eqn.(|l0|) reduces to an identity) @ We follow the conventions by Thorne, Misner and Wheeler a a 3-f = 87rG(p m + Pg ), (12) a' 2 + — r = 4:nGa 2 (p m -p m + p g -p g ), (13) where p g and p g are the averages over nonlinear terms. It is convenient to interprete p g ,p g as the effective density and pressure of the gravitational radiation field. Sub- tracting the averaged field equation ( p~3| ) (multiplied with 5ik) from ( pd| ) and neglecting higher-order terms of the type (A) — A, as is usually done in a self-consistent field approximation, one obtains a modified wave equation for the amplitudes hik ik Ah lk - 2—h lk + 2h ik ( a' 2 a with b = 4TT(Ga 2 (p g ~ Pg ). + 6) = (14) (15) This equation differs from the usual form by a time- dependent term b representing a back reaction of the energy and pressure of the waves on their propagation (a factor of £ was introduced in front of b in order to switch off the back reaction term for comparison pur- poses) . Since p g and p g depend on solutions of the wave equation, some nonlinearity is thus introduced, b re- sults from second-order bilinear terms and is therefore neglected in linear treatments of cosmic wave propaga- tion, but should be kept in the spirit of our approach. In situations where the waves do not appreciably influence the scale factor evolution, the back reaction term can be neglected. Also its influence is usually small in the high-frequency approximation, when the wavelengths are small compared to the Hubble distance. Note also that the back reaction term is exactly zero, if the wave back- ground has the Zeldovich equation of state p p — p g . On the other hand, if low-frequency radiation contributes ap- preciably to the average density and pressure and hence to the scale factor evolution, the modified wave equation must be considered in general. (A previous paper ( |17|) on the same subject was based on the linear wave equa- tion with b = 0. We shall continue its use in Sec. IV in order to compare later the results with those of the general case b ^ 0). Some authors use different definitions of the wave am- plitudes by applying factors of a on them. We have de- fined the spatial components hik as perturbations to the three-metric a 2 Sik- Equivalent to (111) are (h lk /a 2 ) -A(h tk /a 2 ) + 2-(h lk /a 2 ) + 2bh lk /a 2 = a (16) and {hv, /a) -A(h lk /a) + (2b ){h ik /a) = 0, (17) a which are sometimes easier to use. It was fre- quently noted that gravitational wave perturbations in a Friedman-Robertson- Walker (FRW) universe may be described as a pair of massless minimally coupled scalar fields in the same background space-time (see, e.g., Ford and Parker in ( M). If back reaction is important, this correspondence is lost in general, but for a de Sitter scale factor a ~ 1/rj the modified wave equation has the form (□ + £R + m 2 )(<$>/a) = with £ = | - ±m 2 aV, which characterizes non-minimally coupled and possibly mas- sive scalar particles. To calculate the components of the gravitational stress- energy tensor, we assume that the random process rep- resented by hik is homogeneous and isotropic. The cor- relation functions which enter (0,0) are related to cer- tain spectral densities. To see this, we represent hik as stochastic Fourier integral (see B for a detailed treat- ment of the spectral representation of random processes) : hik(x,rf) = j lk (k,r))e lWx dk + conj.compl. (18) From (|l4|), the Fourier amplitudes satisfy the ordinary differential equation 7a " 2 77- fc + 7^(2^ -2^- + k 2 + 26) = 0. (19) Amplitudes of the correlation functions may now be writ- ten as frequency integrals over spectral densities: (^ fc (x, r))hi m (x, r])) = / ((jikim + "f* k i m ))dk, (20) whereQ §} (7ifc(k,r?)7f m (k,r?)) = <5(k - k) 7ifcim , (21) (22) It is easy to see that the spectral densities jikim satisfy the symmetry relations * liklm — flmik i ^iklm — ^kilm j likll = 0, Jiklmk" 1 = 0. (23) (24) The general solution of the algebraic constraints ( 23[ 24) for the spectral densities contains four complex functions of k and rj, which describe polarized background radi- ation in general. The stochastic background of gravi- tational waves expected from pre-galactic stages of the Universe could be polarized due to strong anisotropics expected at Planck time. Here we confine the discussion to unpolarized radiation, represented by a single (real) spectral function a(k,rj). Then "/ikim can be written as: liklm = Ct{k,T])Siklm, Siklm = SuSkm + SimSkl ~ 5i k 8l m , (25) (26) ' We refer to equal space - equal time correlators of the form (poh as "correlation functions". (27) where Oik = uik fcikk j «* 7 & = K K . The transversal <5-symbol satisfies 4 = 2, 4fc fe = 0, 44; =4- (28) To calculate the averaged stress-energy tensor, we need also the spectral densities (likllm) = P s ( k ~ k)5 ik im, {likl'Cn) = 7^(k - k)Siklm, (29) (30) where /3 is real and 7 is complex in general (the correla- tion functions a, fi and 7 were introduced in J14|). Simi- lar to (pO) , the correlation functions which are needed to calculate the effective wave stress-energy tensor can be written as \'lik,r'llm,s) — ^ I ^Oiklm^r^s^'^^ \llik'l'lm,rs) — \i^ik,ri^lm^s) 5 ( h ik h 'lm) = 2 / PSiklmdk, (h ik ti lm ) = 2 / 5R(7)4 im dk, {hikhim.r) — 2 / a/c r 4-; m dk, (hikh' lmr ) = 2 ^(rfSikimkrdk {hikhim.r) — ~\hikhi mr ). Since a, /3, 7 are functions of k (and 77) only, the angular integrations can be performed easily. This gives with Piklm = 44m + SimSkl- (hikhim) = -——{ipikim - 2S ik 6 lm ) / ak 2 dk, (hik,rhlm.s) — 77— ■(— 10&tfe4n4s + HPi/rim4s 105 ^Pilrs^kra ^Plkrs^ira dPkmrsVil ^Pimrs^kl (hikhi. (h'ikh'i m ) = -r--(3pifej m - 24<5; m ) / /3£; 2 dfc, (h lk ti lm ) = -— (3p jfc z m - 24<W / K(7)fc 2 rfA;, (hikhim^) = 0, (hikhlm,r) = 0, {hikh' lmr ) = 0. Since only the real part of 7 enters the averaged ex- pressions, we denote subsequently this real part by 7. Some unusual results can be expected when one deals with stochastic averages over nonlinear quantities. For the averaged three-dimensional Ricci tensor we obtain, using the just derived relations, <^4 3) > = |^4 J ak A dk. (31) Thus, whereas the three-dimensional background met- ric is flat, a random superposition of gravitational radi- ation to this background produces an averaged 3-Ricci tensor with positive curvature, if curvatures are defined as eigenvalues of the Ricci tensor. On the other hand, the stochastic average of the 3-curvature scalar R^ = g lk R ik is negative: (ij(3)) = _??E f a kUk. a b J (32) The bilinear terms which enter R^ and produce a nonzero stochastic average are different from those en- tering Rrj^, thus the result (Ri k )(g %k ) ^ {R lk g lk ) docs not come as a surprise. But it makes it hard to inter- prete stochastic averages geometrically. Energy density and pressure of the gravitational waves are given by similar integrals Pg = 2Ga 6 Pg l-g [ dkk 2 {k 2 a + /3 + 4 7 - - 12cA-), (33) \— [ dkk 2 (7k 2 a-5[3 + 20j—-20a'^-). (34) xa b J a cr 6Ga 6 The deviation from the HF equation of state follows as p g - ip g 1 Ga 6 ,'2 dkk 2 {-ik 2 a + 3/3 - 87— + 4a-5-), and the back reaction function b is given by 8ttC 3a 1 dkk 2 (-k 2 a + 2/3-2~f- 4a- (35) (36) The time evolution of the spectral functions a, /3 and 7 is obtained from ©. Differentiating @,(§9]),(|^) with respect to 77 and using (|19|) leads to a system of coupled differential equations: 27, .'- = 4-/3 + 2 7 (2- - 2\ - 2b - k 2 ), a a a 11 7 a , a 2— 7 + /3 + a(2 — a a 26 -r (37) (38) (39) With given initial values the time evolution of all spec- tral densities and hence correlation functions follows from this system, if the scale factor is known. The system for a, /3, 7 is actually nonlinear due to the presence of the b term. It can be transformed into a single nonlinear dif- ferential equation for a alone. Substituting a'/2 for 7 in the last two equations and introducing a new function e instead of by means of 11 e = /3- a 4a' one obtains , , a a. e'+ 4— )e = 0. a a This is integrated to e (k)c (40) (41) (42) Then the basic spectral function a satisfies the nonlinear differential equation over the spectral density f(k,rf), denoted as "moments'' subsequently: /o = / dk k eo(fc) f(k,r,Y h dk k 2 nM) h = / dk k 2 f(k, V ), h= dk fc 4 /(M). (48) In the density and pressure equation, / and /1 appear only in the combination .94 = /o + A/4. We also note the useful relations ffi = -/i+(--2&)& a /£' = 2.g 4 + 2(— - 26)/ 2 - 2/4 (49) (50) (51) 2aa" - a' 2 - 4e a 4 - 4-aa' - 4a 2 (2 2— - 26 - fc 2 = 0. (43) If we define wave amplitudes with a different power of a, for instance as in (|l7|), this equation can be simplified. With a = fa 2 , VI" - f' 2 + 4/ 2 (k 2 + 26 - — ) - 4e = 0. (44) a This equation is also nonlinear in /, but its nonlinear- ity arises because we deal with the spectral density of an expression which is quadratic in the random wave ampli- tudes hik ■ The solutions of (M) are related to those of a differential equation, which is linear in the wave ampli- tudes. Let hi, hi be a pair of real solutions to such that h" + h{k 2 + 26 - a" /a) = h\h' 2 — h%hi eo, (45) (46) where eo is a time-independent function of k only. Then / = h\ + h\ is a solution of (|44| ) . Not surprisingly, the lin- ear differential equation for h is the Fourier transform of the wave equation ( |l7|) for a single realization of the ran- dom process. We also note that the solutions of (|44| ) can be reduced to solutions fhom{k,rj) of the homogeneous part of ( f44|) (eo = 0): Any solution of the full equation ( [44] ) can be written in terms of a suitable function fhom as / = fhom +£o(k)f ho di] fi home (47) It is convenient to write the expressions for p and p and for other stochastic averages in terms of four frequency- independent, but in general time-dependent, integrals which follow from differentiating f\ once and fi twice and using the differential equation for /. For a general scale factor, energy density and pressure may be rewritten as P 9 = ^(U + 54 + 3-/^-7^/2), (52) zGa 4 a a z 3p ^2ij (7/4 - 5 - 94+5 7 /2 ~ 5 ^ /2) The back reaction function 6 is given by *=It(-/4 + 2.94 + -f 2 -^h\ 6a A a a z (53) (54) thus the relations (pO) , (pll) for the moments become non- linear in general. For completeness, we give the ensemble averages of some other geometrical quantities. The mean value of the four-dimensional Ricci scalar can be written as (R) = K + ^( 3 /4 - 3 34 + ^h + -&\ (55) a° a* a a where the first term is the background contribution. The second term may also be written as 8irG(3p g — p g ). Since R = in vacuum, also the stochastic average of R must vanish. Then ( J55[ ) is equivalent to a combination of (|1J) and (O) (with p m = p m = 0), as one can check easily. There are only two independent and in general nonzero components of the averaged Riemann tensor: (i?oioi) = a' 2 - aa" + ^(g 4 + ^/ 2 - -/^), (56) 3 a z a \Tx a' 2 . a' (Ri2i 2 ) = a' 2 + —(3/ 4 - ff 4 - — h - -fit)- (57) 3 a z a The one independent component of the averaged Weyl tensor (O)10l) - (Cl212) = ^(/4 - ^h + f&)- (58) d a z la contains, of course, no background term. All curvature components result from averaging nonlinear terms, they have no direct connection to the first order perturbations of the curvature tensor, whose average is zero. It is im- mediately seen that {g a ^){C a ij,0v) differs from zero for some components, while, of course, the original expres- sion g a ^C ai _ t p L , and hence its average vanishes. Averaging the four local curvature invariants connected with the Weyl tensor is more complicated: Apart from the cor- relation functions listed above further expressions such h ik (x i ) = h ik (x i )-a 2 (e k + &) (61) as (h ij,kl "mn,rsj must be calculated, which introduces higher-order moments such as fa = Jdk k 6 f(k,rj). We summarize the main results of this section. In or- der to describe the back reaction of gravitational radi- ation propagating in a flat Friedman universe upon the scale factor, we have replaced the Einstein field equations by a system of equations for averaged quantities. This system consists (i) of the Friedman equations (|l^,(|l3|) for a cosmic fluid with density p m and pressure p m , to- gether with the corresponding quantities p g and p g for gravitational radiation given by (|5^) , ( p3| ) , and (ii) of the differential equation (H) for the spectral density f(k,r]). The latter equation may be replaced partially by differ- ential relations for the moments of /. The aim is to solve these equations for the scale factor and for / simultane- ously. One has to keep in mind that the results are only reliable if terms higher than second order in the original field equations can safely be neglected. III. GAUGE TRANSFORMATIONS We have chosen to work in a fixed gauge. As re- cently noted by Unruh llq], there is no basic difference between a gauge invariant approach, where the equa- tions are written in terms of "gauge-independent" vari- ables (see, e.g., [Oj), and an approach where the gauges are fixed in some way. It is only important to fix the gauges completely, leaving no room for residual gauge freedom, which could be mistaken as physical degree of freedom. We therefore ask for the coordinate transfor- mations x^ = x^ + £ M preserving the eight constraints ^oo = 0, hoi = 0, ha = 0, hik.k = 0. To linear order in £ M the perturbations transform as V(»") = M*") - &? $}& ~ fl$&- (59) A simple calculation shows that for a general scale factor the transformations we are looking for must satisfy e° = o, e = e{x k ), e, 0, £ u- 0. (60) The functions £ M therefore do not necessarily form a Killing field, and the first-order perturbations are gauge- dependent in general. The spatial components transform as It is easy to see that the first-order wave equation ( fI4| ) is invariant with respect to (pTI).The transformation of the correlation functions is derived assuming that the vector £* is also a random process with zero mean and indepen- dent of the wave process. We outline the procedure. The random process £* must be homogeneous and isotropic as hik, so its correlation function can be written (ee dk(Fi(fc) k l k h F 2 {k)5 tl (62) £\ = is translated into Fi(k) = —F 2 (k) in the spectral region. It is more difficult to write the condition ^ l kk = directly as condition for the spectral density F\, so we first proceed without this constraint. With ( pl|) one forms (hikhim), this gives J 2 = f 2 + — I dkk i F 2 (k). (63) Similarly, the reduction of (h ik h lm ) yields g 4 = 94 + ^-Jdkk 4 F 2 (k). (64) Writing out (hik, r him jS ) shows that the average of the ^-depending terms does not have the symmetries of the other terms. This reflects the fact that the condition £ l fcfc = was not taken into account. We conclude that the spectral density satisfies the integral relation J dkk 6 F 2 (k) = 0, and that f 4 transforms as h = h (65) It is easily seen that the density p g and pressure p g as defined by (]35 ),(p4] ) are invariant against the transfor- mations (|63])-(35|), and this holds also for the averaged components of the Weyl tensor and for the Ricci scalar. As is well known, the effective stress-energy tensor is not gauge invariant with respect to general coordinate transformations, i.e., transformations which violate the constraints. However, as shown by Abramo, Branden- berger and Mukhanov ( JO] , see also Q , pj), the gauges change at the same time the background geometry to sec- ond order, and these changes just compensate the change of the energy-stress tensor. IV. SOLUTIONS INCLUDING MATTER FIELDS We now try to find solutions of the equations derived in section II, assuming that apart from gravitational ra- diation other forms of matter are present, which domi- nate dynamically during most stages of the cosmic evo- lution, if not always. Dominance of matter means that the 6— term in the wave equation can be neglected. A. The high-frequency regime The high-frequency (HF) regime is defined by the as- sumption that all wavelengths are small compared with the temporary Hubble distance. Then the time deriva- tives which occur in (|37|-p9|) are small: Let A be one of the quantities a, (3, 7, then A' is of the order A/T, where T is a Hubble time, T ss a/a' , and this is much smaller then a multiplication with the wave frequency k: A/T -C Ak. The terms (3 and k 2 a in (B3) are therefore much larger than other terms and must cancel in a HF approximation: /3 = k 2 a. (66) Using this relation and neglecting time derivatives in (|3£ p4|), one obtains Pg = 3Pg = 1 Ga 6 dkk a. (67) equivalent to the equation of state for a gas of noninter- acting massless particles. A similar cancellation of terms must occur in (|43|). This requires £q > 0, gives /eo<2 /k (68) and identifies the function eo(fc) as closely related to the time-independent spectral function in the HF regime. The energy density in this regime, 1 Ga 4 dkk -Jcq, (69) shows the typical 1/a 4 dependence on the scale factor for background radiation. Apart from an overall redshift factor, the gravitational wave spectrum does not change in time. The expression ( p7j ) for the effective gravitational en- ergy density and pressure of high-frequency gravitational radiation holds also in the case that this radiation con- tributes appreciably to the background geometry, beside other forms of matter. The result is only a changed time dependence of the scale factor a. B. Stress tensor evolution in a fixed background Considering now the full spectrum including both low- frequency (LF) and HF modes, one must return to the solution of the complete equation (W4) . The change of the spectrum during expansion depends on the wavelength. The solutions / of (B3) for power law scale factors can be represented by Hankel functions of second kind, if we neglect the b— term assuming that gravitational waves do not contribute appreciably to the scale factor. We prefer explicit expressions, since this allows to use Fourier methods. For the matter universe (a ~ rj 2 ), the radiation universe (a ~ 77) and the de Sitter universe (a ~ 1/rf) as background geometry the general solution can then be written a /a 2 = f = 2npp* + (I + im)p 2 + (I - im)p* 2 . (70) /, m, n are three functions of k and connected with eo by e = 4k 2 (n 2 -l 2 -m 2 ). (71) p(x) is a complex function of x = krj, given by p(x) = (1 + is/x)e ix (72) with s = 1 in the matter and de Sitter universe and s = in the radiation universe. (In spite of the different scale factors, the spectral density / in the de Sitter cosmos has the same form as in the matter universe, since a" /a is the same. This coincidence holds only for the spectral density / and is destroyed, if one returns to a). For large x, f reaches the same asymptotic form in all models. This asymptotic form coincides with the HF expression / = y/eo/fc considered previously, if the ratios — and — tend to zero for large k. How do the different wavelengths contribute to the ef- fective stress-energy tensor? It is usually assumed that modes with wavelengths larger than the temporary Hub- ble radius can be discarded, since they are unobservable and not true waves Efj| , have no appreciable influence on the local wave energy density pl|, |22f|, [J23[ , or appear locally as gauge transformation |24j| , |18||. The relations (p2p and (p3) allow a definite answer. For a radiation universe a ~ 77 one obtains: 7 7 3 a A Gp g = 2n A -n 2 - jj-^Vkfa) + ^Z^M, ^ 5 5 5 3a 4 G Pg = 2n 4 - -^n 2 - l^Mv) + T^^iv) -\^{jl) (74) with the time function ^ 2 {v) = 2 / k 2 (l(k)cos{2kri) ~ m(k)sin(2kr]))dk, (75) Jo we also use the notation ni — L n(k)k l dk for the time- independent moments of n(k). The spectral functions l,m,n are constrained by the condition that the fre- quency integrals (including the moments n 2 and 714) con- verge. To justify the assumption of a pre-determined scale factor, the amplitudes must be small enough that their contribution to the background geometry is negligi- ble. If the function ip 2 (77) with its first time derivative is bounded for large r\ and ip 2 tends to zero for large 77, the relation p g = p g /3 between pressure and density follows asymptotically for large 77, as seen from (|73|) , (f74[) . For small times however, p g and p g show strong deviations from this relation. In particular, p g as well as p g become negative for sufficiently small n. In the case I — m = (and therefore ip2 = 0) this is immediately seen from ( 173 ) and (um . p g and p g are also negative for small n in more general spectra. The main reason is that the spectrum for wavelengths exceeding the local horizon scale at n (corresponding to frequencies k < 1/n) gives a negative contribution to the total values of energy density and pressure. If one goes back in time, the superhorizon modes ultimately domi- nate the spectrum (we always assume that the spectral density in the ultraviolet region drops sufficiently fast to ensure convergence for k — > oo. Ultraviolet divergence is a problem of quantum theory and does not concern us here) . Since the HF part of the spectrum satisfies a p g a 4 = const-law and is therefore not very interesting, we con- centrate on the LF tail, which is essential for back reac- tion effects. C. Back-reaction solutions in the regular low-frequency limit In the long-wavelength or LF limit JnOwe assume that the spectral density / can be expanded in powers of k around k = 0. This assumption excludes a singularity at k = 0, which in some cases can be interpreted as lead- ing to a finite infrared contribution to p and p (Section V). Writing down (44) explicitly requires to calculate the back reaction term b and hence the integrals in (Hq). In the sense of our approximation, this is trivial if the in- tegration is assumed to extend to a maximal (but still small) k = k\. Expanding also eo, (J44]) leads to 2qq"-q n 327rfc|C , 1,2 2 , 1 , 2jL ,a' 6 5 6 a Aqq'- a (76) as lowest-order equation of a hierarchy of equations for the expansion coefficients (we have written q for limfc^o f/a 2 and again eo for lim/j^ofo)- Only this lowest-order approximation will be discussed. Since the solution of ([44 ) can be written in the form (|47|) , it is ap- propriate to treat the case eo = 0. Here q — const is an obvious solution of (|76|), provided we neglect the term proportional to k\ in the bracket, which is small com- pared with the other terms in a LF approximation. The effective energy density and pressure of the extreme LF background follows as Pa 3 _ qk\ (77) Note the a 2 decay compared to the a 4 decay of the energy density of the high-frequency radiation field (cf. |12[ ) . Thus the influence of the LF gravitational radia- tion field on the scale factor is small for sufficiently early times compared to other forms of radiation. The grav- itational wave background considered here is alone not able to support a cosmological model: Since the time or scale factor dependence of energy density and pressure is already fixed, the two equations ([l2]) and ( |l3| ) are not compatible. But we may assume that a relativistic fluid with a p m — p m /3 equation of state is also present. This ensures compatibility at the price of fixing the density ratio between fluid and waves and gives _ 15a' 2 - 4a 2 fcj 5 07r a" + -nkfqa = 0. 5 From the second relation the scale factor follows as a = aisin(lrj), I = 4irk{q/5, (78) (79) (80) thus the energy density of the relativistic fluid changes as fcfg(3-4sin 2 (/?7)) Pm — 10Ga 2 sin 4 (Zr/) (81) It appears that the mere existence of a LF gravitational wave background could cause a reversal of the expansion, however small its energy density may be initially (that is, for times n with In << 1) compared to that of the rela- tivistic fluid, note that p g /p m = sin 2 (In) / (3 — 4 sin 2 (In)) . For In > 0.8861, p g exceeds p m , for nl > it/3 the fluid energy density becomes formally negative through the interaction with the gravitational wave fields. The aver- aged Weyl tensor given by (C0101) — (C1212) — lQi:a\k\q 15 sin 2 (In) (82) vanishes at the singularity n — and is non-zero for n > 0, suggesting that the presence of a LF component of gravitational radiation is not a coordinate effect. The case of a nonzero limit eo for k — > is more complicated and seems not to allow an analytic solution. As another example we treat the interaction of the LF background with a A— term. It is common practice |25| 1 to write the A— constant on the right-hand-side of the field equations, i.e. to consider a fluid with the com- ponents p m = g^g, p m = -gTg. Considering again the case lim fe ^ e o = 0, we arrive at ([77|) as before, but the field equations ( |l2| ) and ( |l3|) lead now to an expo- nential increase or decrease of the scale factor (in the 77-coordinate) : a\ exp (±rjJ Aiik\q/3 (83) As in the former case of a relativistic fluid, the "A— matter" and the LF gravitational wave background are coupled since 45GA (84) y 64 7 r 2 fc 1 13 is required for consistency. The Weyl tensor components I6na 2 klq (C 0101/ (C : 1212/ 15 (85) are also finite (and nonzero) at ?/ = 0. V. DOMINANT GRAVITATIONAL WAVE BACKGROUND Let us now assume a pure gravitational radiation uni- verse. In the high frequency approximation all is already done: We do not have to face an appreciable back re- action term in the wave equation, the spectral shape of the waves is time-independent, the amplitude decreases only due to redshift effects, the equation of state can well be approximated by p g — 3p g , and (|12|) , (|13|) (with Pm = Pm = 0) give the Tolman radiation cosmos. The presence of only LF components in the spectrum does not allow a solution based on gravitational waves alone, thus the spectrum will not be restricted subse- quently. We use the two equations ([30]) and (|5l]) for the frequency integrated quantities g^ and f 4 . Adding (O),(03) with p m — and p m — 0, one has two further equations to determine the four unknown time functions /2,/4 7<?4 and a. We solve the last equations for g^ and /4 and use the result in (pl|) to find / 2 . This leads to g 4 = — (aa" + 3a' 2 ) + — — / 2 - ~-/ 2 , 87r 6 a z 6 a 1 / n n m 1° a ' 2 r /* - 5-(-<w» + 3a' 2 + ^r— h sir 6 a z 5 a' 3 a /2 2- (86) (87) 1 /^^" + (2(1-2 C )- + -(1-6C)^)/: + 3a /2 ' = 2^((3C + 2)^ + (3C-4)^)/ 2 -,/2 + ((C-2)^ + (C + 4)^)/^ a a^ (89) This system may be studied in two cases, (i) assuming C = 0, that is, no back reaction in the wave equation and (ii) full back reaction, C, = 1. Furthermore one has to solve the generalized wave equation ([14]) and check the compatibility of its solution with the moments derived from (|8q)-(|89|). Finally, if all this succeeds, the energy density and pressure of the waves follow from the familiar equations p g = 3a'7(87rGa 4 ), p g = (-2aa" + a' 2 )/(87rGa 4 ). (90) A. C, = 0: Tolman universe With C, — 0, ( JS9[ ) can be written as the product of two factors: (2a'/ 2 - a/a) (aa" - 2a' 2 ) = 0. (91) Thus two different cases emerge, depending on which fac- tor vanishes. If the first factor in (^Tj) is zero, one obtains with an integration constant c (which may be gauged to zero) h = ca 2 (92) This relation is compatible with ( |88| ) if and only if a" = 0, that is, if the scale factor has the time dependence a = brj of the Tolman radiation cosmos J26| , with the energy density and pressure given by p g = 3p 9 8irGb 2 r] 4 (93) From (Bq, pTT) it follows that the other moments are time- independent: .94 = b 2 (c+ —), / 4 3^ 8tt' (94) The average of the Weyl tensor component Cbioi or C1212 comes out as the constant 2b 2 . The expressions for /a, /j and 174 found as solutions of differential equations must be compatible with those derived directly from the the spectral density ( |70[). The three frequency dependent functions l,m,n in ( |70| ) must therefore be chosen so that the moments ji,f^-,g\ have the time dependence which we have just derived. Self-consistency is only ensured if the functions l,m,n exist. Working with real quantities, the spectral density / is given by / = 2(n + I cos(2krj) — m sm(2krj)) (95) in the radiation cosmos. As one verifies from the defi- nition of 54 in jigj ) and of f and /1 in (jig), 34 can be written .94 = 4n 4 - / 4 . (96) From the definitions of / 2 and f 4 in (fig ) one obtains together with (|f) and (fj) k 2 (lcos(2krf) - m sm(2kT]))dk = cb 2 f] 2 /2-n 2 , (97) k 4 (lcos(2kr]) - msm(2kr]))dk = -cb 2 /A. (98) ( pq ) can be obtained from ( |97| ) by differentiation with respect to 77, so only ( B7P is needed. All fc-integrations considered so far run from fc = to infinity. We formally extend l(k) and m(k) to negative values by l(k) = l(-k), m{k) = -m(-k) for k < 0. (99) This allows us to rewrite ( |97| ) as complex Fourier trans- form i '2nJ- (f)(k) exp(2i/c?y) dk 2^2 2 c6 ry 7T 2~ "-2J (100) applied to the complex valued function <fi(k) — k 2 (l(k) + im(k)). Since 4>*{k) = <p(—k), the Fourier transform of 4>(k) is real. To find <f>(k), one has to invoke the Fourier inversion theorem. This requires to consider the right- hand-side of ( |10C| ) also for negative values of rj. Since the right hand side is not absolutely integrable over the whole time axis, <p(k) must be understood as generalized function (see, e.g., J27), p8| , [ g9j for a confirmation of the subsequent calculations). Extending the function space in this way, the Fourier inversion theorem remains valid. For instance, a polynomial in r\ gives rise to the Dirac delta function S(k) and its derivatives in Fourier space. One obtains k 2 l{k) cb 2 d 2 S n 2 4 dk 2 4-7T S(k) k m (k) = 0. (101) (102) We may also extend n(k) to negative values by n(k) = n{—k) for k < 0. The spectral components of energy density p(k,rf) and pressure p(k,rf) are then symmetric functions of k, and the frequency integrated total density can be written p = j/_ dkp(k,rj), with a similar ex- tension of the integration interval for the pressure. This allows to handle terms involving delta functions apply- ing the usual rules for these functions. Replacing second derivatives of the delta function using the formulae (cf. p7[ , the prime here denotes the derivative with respect to k) s(k)8"(k) = s"(0)S(k) - 2s'(0)S'(k) + s(0)S"(k), (103) one obtains for the spectral decomposition of the energy density and pressure from (|52|), (p3l) a 4 Gp 9 (£:,?7) = 2n(fc)fc 4 - -^n{k)k 2 + -^-cb 2 5"{k) '/ (il n2 _ cb 2 )5(k), 7f (104) a Gp g — 3a Gp g 2tlA b^c ~~2 (106) in agreement with (|96|), showing the self-consistency of the calculation. The singularities of the spectral decomposition p g (k, rf) show that the infrared mode k — contributes a finite and time-dependent amount Ga A pi r = -^n 2 — ^|- to the total energy density. At superhorizon scales, more pre- cisely at scales with kr\ < %/3.5, the spectral components become negative. Similar conclusions hold for the pres- sure, which becomes negative for kr\ < V2.5. The in- tegrated values of energy density and pressure however stay always positive thanks to the infrared behaviour of their spectral components. B. ( = 0: de Sitter scale factor The vanishing of the second factor in (pll) gives the scale factor of the de Sitter universe, a = -4-. The density and pressure calculated from (90) are time- independent: Pg = ~Pg = 3H 2 8^G' (107) The differential equation for f 2l (pq), has the general solution t r . 7/3 3 ln M h = -, + soil 1 - r/* UttH 2 if (108) with two constants r$ and sq, where tq may be gauged to zero. From (J8qJ87]) one obtains , _ 65 1/3 27 1 •' 4 = ~7T Sor l ~ 7771 tj2 4 - 9 104 7T_H z! ?7 4 (109) .94 1° 4. ^W/3 _ ±J^± 4- il^_ (HO) '] -soV 13itH 2 t] 4 104 nH 2 ri 4 and the Weyl tensor components arc found as 728 10 1 (CoiOl) = (C iai2 ) = -^ W /3 ~ ^^4- (HI) We have again to check the compatibility of the spectral density / (defined by (f70|)) with the time dependence of its moments f 2 , f±, g^ as given by the last three equations. The integral f 2 corresponding to the definition in (Uq) is 3a 4 Gp„(fc,77) = 2n(fc)fc 4 - %n{k)k 2 + --\cb 2 6"(k) rj z ir] z + {^n 2 -cb 2 )5(k). (105) rj z The spectral quantities can be integrated immediately to give the finite total values f 2 = 2n 2 + -^n + 2 k z (l cos(2krj) - m sin(2kr])) dk V 2 Jo V Jo (—1 cos(2krj) + m s'm(2kr))) dk k(m cos(2krj) + I sin(2kr))) dk. (112) 10 A similar expression holds for fa . Straightforward calcu- lation shows that the integrand of 54 = L k 4 g(k, 77) dk defined by (f49|), has the form (x = krf) g = u cos(2x) + v sin(2.T) + w, (H3) 2 u = —r(-l(x 4 - 3.t 2 + 1) + 2mx(x 2 - 1)), a; 4 v = -r(m(x i - 3x 2 + 1) + 2lx(x 2 - 1)), x 4 w = 2n(l - l/x 2 + 1/x 4 ). Rewriting the just derived integral expressions for fa and #4 as well as for fa as complex Fourier transforms, we obtain: 2 1 2i fa = 2n 2 + -^n a + ijj 2 - -ttVo + —ipi, (114) 77 z ry z 7j 2 1 2i / 4 = 2n 4 + — n 2 + ip4 5-^2 H "03, (115) rj z i~i z 77 2 2 2i , .94 = 2n 4 5-712 + ^-n - ^4 03 77 z r) 77 3 2i 1 H — jV>2 + -o'Vi 700, (116) 7j z ry 3 77* where ipj is a family of time functions defined by (117) ipj = k 3 (I + im) exp(2ikrj) dk (ip2 was already introduced in Section IV) . As in the Tol- man case we have extended the functions I, m, n to nega- tive values of k by l(—k) = l(k), m(—k) = — m(— k) and n(—k) = n(k), thus the functions ipj with even (odd) j are real (pure imaginary) . The members of a ^-family are connected by differentiation and integration according to the rule %//j = 2iip 3+1 , (118) which holds also if the ipj are generalized functions. The aim is again to obtain the complex spectral den- sity l(k) + im(k) from any of the time functions ipj by invertin g the corresponding Fourier integral. We start solving (115) for ip 2 , where fa is replaced by the expres- sion (109). Using the rule ( 118 ) repeatedly, one obtains the differential equation r, 2 iP% - 4777^ + 402 8n 2 + 8n 4 i] 2 + 27 260 267rilV -soV 7/3 (119) for ip2- Its solution is given by (note that adding a solu- tion of the homogeneous equation would give the wrong time dependence) 1P2 = 2n 2 - 477 4 7j 2 + 13s 7j 7/3 + 52ttH 2 V (120) If one member of a 7/;- family is known, other can be found by differentiation and integration. Straightforward cal- culation shows, that the T^-functions derived from (12C) satisfy also (114) and (116), if fa and 54 on the left-hand- sides are substituted from ( |108| ) and (110). We now apply the Fourier inversion theorem to (117) with j = 2 and ip 2 taken from (12C). This requires a continuation of the (real) function ip 2 into the region 77 < 0. Here only the term proportional to sq requires more consideration, but this term contributes nothing to p g and p g , we can there- fore put the integration constant sq equal to zero. One then obtains k l(k) = 77,4 k 2 m{k) = 0. d 2 5(k) dk 2 2n 2 S(k) - m 26ttH 2 ' (121) (122) Again infrared modes enter the spectral density l(k), but contrary to the Tolman case one is free to specify the spectral function n(k). If n(k) is chosen as zero, the delta function terms in / and in the spectral decomposition of p and p vanish, but / has still singular terms seen in an expansion around k = 0: / 1 27/ 4 13nH 2y k ' rfk*' 39nH 2 r<7 + k 3 + o(k 5 ) (123) (we missed in section IVC the de Sitter case, since we had excluded infrared singularities). In spite of this sin- gularity, density and pressure of the gravitational waves integrate to the finite constant values ( |107| ). Additional terms from a spectral function n(k) 7^ add nothing to the total values p g and p gi thus the singular infrared (k = 0) component in p g is cancelled by the integrated contribution of 77-modes with k ^ 0, the same holds for the pressure. It is easy to extend the calculation to account for a cosmological constant by adding matter terms p m — A/(87rG) and p m = —p m to the equations ([l^Jl3] ). Re- peating the calculation at the beginning of this section, one obtains (for Q — 0) the same product (|9l]) as in the absence of a A constant. The vanishing of the first fac- tor gives a Tolman-de Sitter model, where gravitational radiation has an equation of state (EOS) of the form Pg = Pg/3 an d an energy density decaying as a -4 . The scale factor 0,(77) can be expressed in terms of an elliptical integral. More interesting is the model corresponding to the van- ishing second factor, since here the de Sitter scale factor follows. No new calculation is needed: In equation (107) we have to substitute for p g and p g the total values p g +p\ and p g + p\. In the calculations following this equation we must only replace the de Sitter Hubble constant H in all equations by H = H/{\ - 64ir 2 G 2 A/(3H 2 )) 1 / 2 . The de Sitter expansion is generated by two independent "matter" sources, by a genuine A constant as well as by a suitable spectrum of gravitational waves. We may have 11 an arbitrary (but time-independent) mixture of both in- gredients. For sufficiently small A gravitons dominate. If A reaches the threshold A* = H 2 /(64tt 2 G 2 ), p g be- comes zero and turns to negative values for still larger A, to ensure the same total energy density for a different composition. We have not discussed in this article the origin of the primordial wave spectrum, but it should be noted that the expressions for energy density and pressure intro- duced here may be of interest for concrete models. In the case of a quantum origin due to vacuum fluctuations in a de Sitter cosmos Q, the produced gravitons can be described by a two-point correlation function, which for Bunch-Davies vacuum |3l]] corresponds to a classical spectral density / given by /, BD —a 1 kir 2 rfk 2 )■ (124) Ibd is obtained by comparing the quantum expectation values of bilinear terms in the metric (as given, e.g., in |p2fl ) with the stochastic averages discussed here. The spectral densit y Jbd has the same infrared singularity as the expression (|l23| ) derived previously. Since the coinci- dence of both spectra holds only approximately for small k, the back reaction of the Bunch-Davies gravitons on the scale factor will change the scale factor^]. Again, in spite of the singularity of the spectral density /bd, the inte- grated values of energy density and pressure require no infrared cut-off. If one introduces an ultraviolet cut-off at the frequency k\ with x± = kirj, one obtains for energy density and pressure of the Bunch-Davies gravitons Pa HH 4 x\(x\ - 7), p g hH 4 xfixf + l). (125) Considering limiting cases, for high frequencies ii > 1 follows the expected EOS p = p/3, for low frequencies iiCl one has p = —p/i (together with p < 0). The expressions (125) for p g and p g apply only for a de Sitter scale factor. As a result of to back reaction, the local values of p g and p g change with the background geome- try (thus they are, in a sense, not local): The lower the frequency, the stronger is the dependence of the effective EOS on the background gravitational field. C. f = 1: Full back reaction models Unfortunately, the approach described in the previous subsections does not work in the real case £ = 1. Even ^Similar conclusions follow from the important work by Tsamis and Woodard, who have discussed quantum gravity back reaction on an inflationary expansion rate in numerous papers 0. the first step, finding analytic solutions of (|88|) , (|8£|) for a, /2, has been so far unsuccessful. The back reaction function can be written b = C(ir + ^~) ^ or a domi- nant gravitational wave background, thus the differential equation for the wave amplitudes h takes for C = 1 the unusual form n" n' 2 h" + h(k 2 + — + 2%) = Q. (126) This equation suggests that graviton creation persists in a wave dominated universe, provided the scale factor is different from a ~ 7/ 1 / 3 . Note that the latter scale factor is inconsistent with the system (|88|) , (p9|) . For a treat- ment of the general case one has to resort to numerical calculations, which will be discussed elsewhere. VI. FINAL REMARKS We have seen that the stochastic back reaction equa- tions treated here form, on the one hand, an apparently self-consistent system of equations with interesting solu- tions. On the other hand, it is not clear, how far the solutions deviate from true solutions of Einstein's field equations, either for some range of parameters or in some regions of space-time. We shortly discuss what could be done to clarify and to improve the situation. One has to realize that ensemble averages are consid- ered, which always differ from actual realizations of a ran- dom process. This is inherent to the method and cannot be changed, but one is able to say something more about statistical deviations from true solutions, if a Gaussian or some other process is assumed. The main shortcoming of the approach is the use of a second-order approximation to general relativity, but improvements are possible. In principle, Monin and Ya- glom's statistical treatment of nonlinear field theories works for arbitrary nonlinearities, if they are present in polynomial form. Only the technical complexity grows in a full treatment, since many higher-order correlation functions must be taken into account. Writing the Ein- stein field equations as 1 (R^ ~ ^uR pa 9 P <r)g a = KT^g* (127) where g is the determinant of the metric tensor, we see that the left-hand sides consist of huge polynomials of maximal degree 12 in the (exclusively) covariant compo- nents of the metric tensor and its first and second deriva- tives. Taking the ensemble average of these expressions leads to correlation functions up to sixth order. If the random process g^ v is Gaussian, the higher order corre- lation functions can be reduced to the second-order func- tions studied in this article. This would allow us to turn the back reaction equations into -in some sense - exact 12 relations. Present discussions on a primordial stochastic gravitational wave background usually assume a quan- tum origin of gravitons, which are born out of zero-point vacuum fluctuations. The classical correlations discussed here are expected to be related to quantum mechanical expectation values, thus it seems natural to assume Gaus- sianity. P. Anninos, J. Masso, E. Seidel, W.-M. Suen and M. To- bias, Phys. Rev. D 56 842(1997), P. Papadopoulos, E. Seidel, L. Wild, Phys. Rev. D 58 (1998) 084002 [2] [3] [4] [5] [6] ['] [9] [10] [11] [12] [13] [14] C. Bona, J. Masso, E. Seidel and P. Walker, jt- qc/9804052 R.D. Brill and J.B. Hartle, Phys. Rev. 135 B271 (1964). R.A. Isaacson, Phys. Rev. 166, 1263; 1272 (1968), Y. Choquet-Bruhat, Coram. Math. Phys. 12 16 (1969), M.A.H. MacCallum and A.H. Taub, Coram. Math. Phys. 30 153 (1973), P. R. Anderson and D.R. Brill, Phys. Rev. D 56, 4824 (1997). M. Carfora and A. Marzuoli, Phys. Rev. Lett. 53 2445 (1984), S. Bildhauer and T. Futamase, Gen. Rel. Grav. 23 1251 (1991), R. M. Zalaletdinov, Gen. Rel. Grav. 24 1015 (1992), 25 673 (1993), T. Buchert and J. Ehlers, Astron. Astrophys.,320 1 (1997), J. P. Boersma, Phys. Rev. D 57 810 (1998). L. Mandel and E. Wolf, Optical coherence and quantum optics, (Cambridge University Press, Cambridge UK, 1995), M. Born and E. Wolf, Principles of Optics, 6th ed., (Pergamon, Oxford, 1980). A.S. Monin and A.M. Yaglom, Statistical Fluid Mechan- ics: Mechanics of Turbulence, Vol. 2 (The MIT Press, Cambridge, Massachusetts, 1975). N.C. Tsamis and R.P. Woodard, Nucl. Phys. B 474 235 (1996), Ann. Phys. 253 1 (1997), Phys. Rev. D 54 2621 (1996), Ann. Phys. 267 145 (1998), Phys. Rev. D 57 4826 (1998), Phys. Lett. B 426 21 (1998) I. Iliopoulos, T.N. Tomaras, N.C. Tsamis and R.P. Woodard, Nucl. Phys. B 534 419 (1998), L.R.W. Abramo, N.C. Tsamis and R.P. Woodard, Fortsch. Phys. 47, 389 (1999) , L.R.W. Abramo and R.P. Woodard, hstro-ph 9811430| , 9811431. L.H. Ford and L. Parker Phys. Rev. D 16 245, 1601 (1977), B.L. Hu and L. Parker Phys. Lett 63A 217 (1977). E. Flanagan and R.M. Wald, Phys. Rev. D 54 6233 (1996). B.L. Hu, gr-qc/9902064 and references therein. L.R.W. Abramo, R.H. Brandenberger and V.F. Mukhanov, Phys. Rev. D 56, 3248 (1997), Phys. Rev. Lett. 78, 1624 (1997). L.R.W. Abram o, Ph.D. thesis (Brown University, 1997), gr-qc/9709049| . [is; [16. [17; [is; [19 [20 [21 [22 [23 [24 [25 [26 [27 [28 [29 [30 [31 [32 C.W. Misner, K.S. Thorne and J. A. Wheeler, Gravitation (Freeman, S. Francisco, 1973). F.I. Cooperstock, V. Faraoni and G.P. Perry, Mod. Phys. Lett. A 10, 359 (1995), Int. J. Mod. Phys. D 5 375 (1996). G. Dautcourt , Proc. Marcel Grossmann Symposium, to be publishe d, |gr-qc/9712042] W. Unruh, |astro-ph/9802323 . M. Bruni,, S. Matarrese, S. Mollerach and S. Sonego, Class. Quan t Grav. 14 2585 (1997), S. Matarrese, S. Mollerach, |astro-ph/9705168 . Ya.B. Zeldovich and I.D. Novikov, Astron. Zh. 46 960 (1970). B. Allen, Phys. Rev. D 37 2077 (1988). D.M. Tavares and M.R. de Garcia Maia, Phys. Rev. D 57 1 (1998). M.R. de Garcia Maia, J.C. Carvalho and J.S. Alcaniz, Phys. Rev. D 56 6351 (1997). V. Sahni, Phys. Rev. D 42 453 (1990). J.M. Overduin, F.I. Cooperstock, Phys. Rev. D 58 (1998) 043506 and references given there. R.C. Tolman, Phys. Rev. 38 1758 (1931). I.N. Sneddon, The Use of Integral Transforms (McGraw- Hill Book Company, New York, 1972). A.H. Zemanian, Distribution Theory and Transform Analysis, (McGraw-Hill Book Company, New York, 1965). G.B. Folland, Fourier Analysis and its Applications, (Brooks/Cole Publishing Company, Pacific Grove, 1992). L. Parker, Phys. Rev. Lett 28 705 (1972), Phys. Rev. 183 1057 (1969), L. Grishchuk, Zh. Eksp. Teor. Fiz. 67 825 (1974) [Sov. Phys. JETP 40 409 (1975) ]. T.S. Bunch and P.C.W. Davies, Proc. R. Soc. A 360 117 (1978). B. Allen and S. Koranda, Phys. Rev. D 50 3713 (1994). G. Dautcourt, Monthly Notices RAS 144, 255 (1969), Appendix D. G. Dautcourt, Astronom. Nachr. 295 243 (1974). 13